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摘 要 : 本 文 在 激光 测 高 方法 获取 的 100m 233889 € H DEM 数据 基础 上 ,选取 1407 个 不 同 
态 的 撞击 坑 ， 提 取 多 视角 刻画 撞击 坑 形态 的 地 形 因 子 指标 ， 并 进行 撞击 坑 形 态 特征 分 析 ; 在 
此 基础 上 ， 采用 数理 统计 与 地 理 建 模 的 方法 ,重点 探讨 地 形 因子 间 的 关联 性 ,进一步 构建 基于 
多 因子 的 量化 模型 ， 以 解释 地 形 因子 信息 与 撞击 坑 形 态 发 育 的 内 在 联系 。 结果 表明 , 深 径 比 (Y ) 
与 直径 ( X1 ) 、 体 积 ( X> ) 、 圆 度 ( Xs ) 、 体 态 比 ( X4 ) 以 及 坑 壁 坡度 (Xs ) 能 较 好 地 刻画 挤 
击 坑 形 态 发 育 特征 ， 通 过 多 元 线性 回归 拟 合 构建 指标 间 的 函数 关系 ， 且 拟 合 模型 在 0.05 置信 水 
平 下 的 显著 性 F=492.037>>F ， 拟 合 模 型 的 的 显著 性 极 好 
Y = (7.941B -007)x x, + (8.063E - 016) x x, +0.082xx3+0.012xx4 0.007 x xz -0.042 为 其 模拟 画 
数 关系 式 ， 拟 合 优 度 达 0.637。 该 模拟 方法 在 月 表 撞击 坑 形 态 发 育 研究 中 有 重要 意义 ， 是 撞击 坑 
形态 特征 的 数字 地 形 模拟 与 信息 挖掘 的 一 次 有 益 实践 。 


关键 词 : 撞击 坑 形 态 ; 地 形 因 子 指标 ; 撞击 坑 发 育 ; 量化 模拟 
中 图 分 类 号 : P184. 8 
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f 坑 不 仅 是 月 球 表面 最 典型 的 形 貌 构造 单元 ， 也 是 月 球形 态 特 征 演化 和 发 育 研 究 的 直接 切入 


点 。 因 月 表 几 乎 没有 大 气 层 保护 ， 且 无 风 、 水 的 侵蚀 和 地 质 活动 微弱 ， 所 以 撞击 坑 得 以 较 完整 的 


保存 站。 撞击 坑 的 形成 过 程 、 形 态 特征 及 空间 分 布 等 从 不 同方 面 为 研究 月 球 演化 历史 提供 了 线索 和 
方法 MM。 撞击 坑 坑 层 是 撞击 坑 形 态 特 征 表达 的 最 直观 载体 ， 是 月 表 撞 击 坑 形态 模拟 重建 的 基本 单元 。 
通过 对 撞击 坑 坑 展 形态 特征 及 其 量化 指标 关系 的 研究 ， 不 仅 可 以 还 原 撞击 坑 的 形成 过 程 ， 还 可 以 进 


一 步 推 知 月 球 地 质 特 征 、 物 质 组 成 ， 揭 示 演 化 过 程 等 由 。 构 建 月 表 撞 击 坑 形 态 量化 指标 因子 模型 ， 
对 于 月 表 接 击 坑 形态 表达 和 发 育 


探索 不 同 指标 因子 间 的 关联 特征 ， 实 现 指标 因子 的 量化 关系 横扫 
演化 有 着 积极 的 理论 意义 和 实践 意义 叫 。 


` 


id 


将 撞击 坑 深 度 直径 比 作为 撞击 坑 相 对 年 龄 的 判定 依据 ， 撞 击 坑 可 以 分 为 新 鲜 、 年 轻 、 成 熟 和 老 
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年 四 类 四。 郑 丽 基于 行星 撞击 坑 分 形 及 分 维 数理 论 ， 对 月 表 撞击 坑 数 目 与 直径 之 间 的 关系 进行 了 深 
入 研究， 揭示 了 撞击 坑 月 表 地 质 年 龄 复杂 的 量化 本 质 ， 实 现 地 貌 形态 量化 研究 从 简单 的 撞击 坑 直 径 
向 分 维 数 的 转变 1。 严 艳 梓 研 究 了 雨 海 地 区 粗糙 度 及 与 岩 性 的 关系 ， 通 过 Hurst 指 数 的 计算 ， 分 析 了 
月 表 粗 烽 度 分 布 特征 及 其 与 岩 性 和 地 质 年 龄 的 关系 ， 进一步 探索 了 火山 作用 对 月 表 粗 糙 度 的 影响 外。 
王 娇 、 程 维 明 等 基于 撞击 坑 数理 统计 和 比较 分 析 ， 从 撞击 坑 的 位 置 、 大 小 、 形 状 、 坡 度 、 方 向、 中 
央 峰 、 辐 射 纹 等 7 个 视角 52 个 小 类 指标 ,构建 了 月 表 撞 击 坑 全 关系 型 数据 库 ， 实 现 对 撞击 坑 数 据 的 存 
储 、 检 索 、 处 理 和 应 用 外 。 王 心 源 、 吉 玮 等 选择 坑 层 等 效 半径 、 坑 底 等 效 半径 和 撞击 坑 的 深度 3 个 参 
数 进行 了 统计 分 析 和 聚 类 分 析 ， 以 此 来 研究 撞击 坑 形 貌 与 3 个 参数 复杂 程度 的 关系 IM。 陈 曼 等 基于 月 
表 撞 击 坑 的 形态 学 特征 构建 了 描述 指标 体系 ， 对 其 分 别 进行 地 理学 意义 的 描述 和 计算 方式 的 表达 ， 
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在 月 表 形 貌 学 和 数字 地 形 分 析 中 ， 地 形 因子 的 表达 虽 已 经 从 定性 描述 转变 到 定量 表达 ， 但 传统 
的 地 形 因子 和 特定 的 描述 指标 受 限于 提取 过 程 和 方式 ， 仪 对 月 表 撞击 坑 形 态 进 行 综 合 概括 和 单一 刻 
画 522。 此 外 ， 基 于 数字 地 形 信息 挖掘 的 多 因子 联合 定量 表达 以 及 量化 关系 研究 ， 是 反映 形态 特征 
与 发 育 演化 和 揭示 多 因子 与 自然 系统 之 间 规 律 的 重要 内 容 。 

本 文 以 月 表 撞 击 坑 为 基本 单元 ， 以 LORC 的 100m 高 分 辨 率 DEM 为 数据 基础 ， 在 全 月 选择 1407 
个 不 同类 型 撞击 坑 ， 它 们 能 和 科学、 全面、 准确 地 刻画 撞击 坑 形 态 特 征 与 地 形 信 息 ， 探 索 描 述 撞击 坑 
形态 特征 指标 因子 间 的 关联 性 和 表征 意义 ， 构 建 基于 多 因子 的 量化 关系 模型 ， 以 揭示 地 形 因子 信息 
与 撞击 坑 形 貌 发 育 的 内 在 联系 。 


1 ”数据 来 源 与 研究 方法 
1.1 研究 区 及 实验 样 区 


撞击 坑 的 形成 及 发 育 过 程 受到 诸多 不 确定 因素 的 影响 ， 如 陨石 撞击 方向 、 力 度 、 月 表 岩 性 与 月 
二 厚度 等 ， 因 此 在 撞击 坑 形 态 指 标 因子 的 计算 中 ,撞击 坑 的 选择 尤为 重要 ， 对 于 撞击 坑 类 型 的 划分 ， 
已 总 结 出 几 种 学 界 较为 权威 的 分 类 方案 ， 本 文 以 陈 巡 等 对 月 表 接 击 坑 类 型 的 划分 作为 研究 基础 "1。 
基于 类 型 全 面 、 发 育 典 型 等 原则 选择 坑 体 结构 较为 完整 的 撞击 坑 作 为 研究 数据 ， 如 图 1 所 示 。 
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图 1 选择 撞击 坑 在 全 月 球 分 布 情况 


Figurel The distribution of the slected craters on the lunar surface 
基于 现 有 资料 ， 本 文 依据 数据 完整 性 、 实 验 科学 性 和 数据 可 获取 性 等 原则 ， 结 合 撞 击 坑 面积 
小 、 接 击 坑 间 的 空间 位 置 关 系 、 撞 击 坑 的 形态 发 育 特征 等 考虑 ， 在 月 表 选 择 1407 个 撞击 坑 ， 其 中 包 
括 离散 型 撞击 坑 、 群 罕 型 坑 连 坑 、 群 从 型 坑 套 坑 以 及 边缘 退化 型 撞击 坑 四 大 类 型 。 具 体 类 型 如 表 1 
所 示 。 


dd 撞击 坑 不 同类 型 
Table 1 Different types of impact craters 


类 型 山体 阴影 图 类 型 山体 阴影 图 
新 鲜 撞击 退化 撞击 坑 
Ut (the (the 
fresh dispersed 
dispersed craters of 
craters) degradation) 
新 鲜 坑 连 退化 坑 连 坑 
Ut (the (the 
fresh connected 
connected craters of 
craters) degradation) 
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新 鲜 坑 套 
Ji (the 
fresh 
contained 


craters) 


2.2 数据 来 源 


月 球 勘 测 轨道 飞行 器 (Lunar Reconnaissance Orbiter, LRO) 获取 的 100 m 分 辩 率 数字 
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退化 坑 套 坑 
(the 
contained 
craters of 


degradation) 


高 程 模 型 


(Digital Elevation Model, DEM) 数据 。 基 于 基本 地 形 因子 和 撞击 坑 形 态 特 征 实现 撞击 坑 地 形 结 构 的 


划分 ， 作 为 研究 的 数据 基础 。 


2. 3 实验 流程 


机 理 ， 有 效 地 实现 撞击 坑 地 形 因 子 的 遂 选 。 
然后 对 1407 个 撞击 坑 地 形 因 子 和 指标 进行 提取 并 量化 统计 。 
在 此 基础 上 ， 利 用 数理 统计 与 地 理 建 模 的 方法 ， 重 点 探讨 撞击 坑 地 形 因 子 间 的 关联 性 ， 进 一 步 


实现 基于 多 因子 的 量化 关系 模拟 。 


2 ”结果 与 分 析 


2.1 指标 因子 的 选择 


对 于 月 表 撞 击 坑 的 发 育 而 言 ， 


出 能 够 反映 撞击 坑 形 态 特 征 与 发 育 


首先 通过 对 已 建立 的 撞击 坑 形 态 指标 体系 进行 总 结 ， 结 合 撞击 坑 形 貌 形态 特征 表达 和 发 育 演化 


坑 层 是 撞击 坑 形 态 特 征 与 演化 过 程 的 基本 表征 ， 反 映 月 表 接 击 坑 
退化 发 育 的 演变 过 程 。 因 此 ， 本 文 基于 已 建立 的 撞击 坑 形 态 指标 体系 ""， 从 撞击 坑 尺 寸 规模 、 形 态 
发 和 特征 统计 等 方面 对 指标 进行 剖析 ， 兼 顾 其 所 毕 含 的 地 学 意义 和 指标 依赖 关系 ， 科 学 合理 地 遂 选 
特征 的 指标 因子 ， 用 作 构 建 量化 模型 的 基础 。 选 取 的 指标 因子 见 
表 2， 这 些 指 标 从 坡度 、 高 程 以 及 二 三 维 空间 形态 等 多 重视 角 出 发 ， 重 点 刻画 撞击 坑 的 基本 形态 ， 


强 含 丰富 的 接 击 坑 形 态 发 育 特征 信息 ， 接 击 坑 指标 获取 过 程 见 图 2. 


表 2 月 表 撞 击 坑 形 态 指标 


Table2 Indices for morphological of the lunar crater 


HER 定义 计算 方法 指标 意义 图 例 说 明 
直径 撞击 坑 坑 层 直 径 Da = V4Area/r 反映 撞击 坑 最 终 状 态 下 的 直径 图 2a 
ial HE 撞击 坑 接 近 圆 的 程度 C = 4nA/P^ 反映 撞击 坑 边 缘 的 复杂 性 
童 击 坑 坑 层 最 高 点 月 坑 底 最 er 图 2b 
深度 m HH x onn 反映 撞击 坑 最 后 状态 下 的 深度 
体积 aon V = A* R* 0.5 - SUM * cellsize 反映 撞击 坑 体积 的 大 小 


体积 
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撞击 坑 坑 层 最 小 外 界 矩 形 的 
体态 比 R, = WL 反映 了 撞击 坑 的 狭窄 程度 图 2c 
宽度 与 长 度 之 比 

深 径 比 撞击 坑 深 度 与 直径 之 比 H/D, 反映 撞击 坑 的 发 育 程度 
gines 撞击 坑 坑 层 最 大 高 程 与 最 低 反映 了 撞击 坑 在 月 表 垂 直方 向 图 2b 
SUS LIK IZ R=H -H 

EX LH KCMAX KCMIN 的 起 伏 程度 
坑 壁 平均 坡度 撞击 坑内 壁 的 平均 坡度 ArcGIS Slope 反映 撞击 坑内 壁 的 陡峭 程度 


TE; Area 表示 撞击 坑 的 坑口 面积 ，P 表 示 撞 击 坑 坑 唇 周 长 ，Hkcmax 和 Hkdmin 分 别 表示 撞击 坑 坑 唇 


最 高 点 和 撞击 坑 坑 底 最 低 点 ，SZ 为 坑 壁 和 坑 底 所 向 下 对 应 的 体积 ，ce77size 表示 DEM 数据 的 空间 
DIR, MERERI RRE WA E, L 表示 坑 唇 最 小 外 接 和 矩形 的 长 。 


Note:Area is the area of the crater pithead, P is the perimeter of the impact crater , Hkcmax and Hkdmin 


respectively indicate the highest point of the impact crater rim and the lowest point of the crater bottom, 
and SUM bounded by the crater bottom and walls, cellsize is the spatial resolution of the DEM data, W is 


the width of the minimum bounding rectangle of the impact crater rim, and L is the length of the minimum 


bounding rectangle of the impact crater rim. 


i b 


图 2 撞击 坑 形 态 参数 指标 的 获取 过 程 


Figure 2 Acquisition process of crater morphological indices 


YE: (a) RAN ALE (Da) 的 获取 过 程 ，(b) 是 深度 与 坑 层 起 伏 度 的 获取 过 程 ，(c) 是 体态 比 的 获取 过 程 。 


Note: a is the acquisition process of diameter, b is the acquisition process of depth and rim fluctuation, c is 


the acquisition process of posture ratio. 


2.2 指标 因子 的 量化 计算 


在 研究 中 ， 首 先 从 刻画 月 表 撞 击 坑 特 征 和 自然 表达 效果 角度 对 现 有 撞击 坑 


BRAKE 


行 对 比 。 结 果 发 现 ， 基 于 地 形 分 析 方 法 的 撞击 坑 坑 层 的 提取 算法 ， 在 一 定 尺 


A 
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取 算 法 的 本 质 认识 发 现 ， 目 前 大 多 接 击 坑 坑 层 形态 的 提取 算法 对 其 识别 准确 
RE 且 伴 随 有 因 窗 口 平 滑 而 造成 的 细节 简化 问题 ， 造 成 撞击 坑 真 实 坑 层 信 
识别 的 撞击 坑 坑 层 形态 :的 提取 方式 对 于 撞击 坑 特 征 的 反映 具有 较 高 的 可 对 ETE, 
D A ASUR RAS A E. dk 1407 个 撞击 坑 样 本 的 基础 上 ， 利 用 DEM 数据 ， 
行 撞击 坑 坑 层 形态 判 识 ， 以 此 来 提取 撞击 坑 坑 唇 线 。 然 后 根据 撞击 坑 类 型 划 
个 地 形 指标 因子 进行 量化 计算 ， 各 项 指标 计算 结果 见 表 3。 
表 3 不 同类 型 撞击 坑 形 态 指 标 平 均值 统计 对 比 


度 都 受 实验 样 区 选择 的 
息 的 丢失 ， 而 基于 目 视 
提取 结果 与 月 表 撞 击 
采用 目 视 识别 的 方法 进 
分 结果 ， 对 已 选择 的 8 


Table3 Comparisons of mean values of morphological indices of craters of different types 


类 型 ”直径 (Dykm) ”深度 (H/km) 体积 (V/km) WE (C) 体态 比 (Rp) — 深 径 比 (HVDJ) SEE CSF) 
1 12.396 2.313 167.554 0.9728 1.1564 0.1867 21.40 
2 11.395 1.919 161.681 0.9356 1.0274 0.1714 17.34 
3 31.199 3.346 2710.198 0.9522 1.1289 0.1122 17.28 
4 47.936 2.135 1226.035 0.9455 0.8838 0.0771 13.48 
5 29.659 2.721 1937.222 0.8980 0.8901 0.1136 14.09 
6 59.432 3.921 8366.350 0.9092 0.8771 0.0750 13.57 


TE: 1 表示 新 鲜 离散 型 撞击 坑 ，2 表示 新 鲜 坑 连 坑 ，3 表示 新 鲜 坑 套 坑 ，4 表示 退化 离散 型 撞 了 
表示 退化 坑 连 坑 ，6 表示 退化 坑 套 坑 。 


Note:1 is the fresh dispersed craters; 2 is the fresh connected craters;3 is the fresh contained craters; 4 is the 


Hbi, 5 


dispersed craters of degradation;5 is the connected craters of degradation;6 is the contained craters of 


degradation. 


从 表 3 WDE HA RIA EH 


fi 坑 态 指标 之 间 的 差异 性 ， 新 鲜 撞击 坑 的 圆 


说 明 新 鲜 撞击 坑 的 坑口 跟 接 近 圆 形 ， 在 新 鲜 撞击 坑 中 ， 坑 连 坑 


态 由 于 相互 的 挤 压 发 4 


间 存 在 相互 


影响 ， 


E 本 本 


[xr 


E 较 大 的 变化 ， 新 鲜 撞 


fieri 


2.3 指标 关联 性 分 析 
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TEASE ME ze 
异性 ， 但 各 要 素 之 间 是 普遍 存在 联系 的 而 绝 非 孤 立 的 
并 实现 撞击 坑 形 态 的 发 育 模 拟 ， 是 本 研究 的 
关联 性 分 析 也 即 相关 性 分 析 ， 是 揭示 描述 地 3 
R。 本 文通 过 所 选择 的 撞 了 


圆 度 


是 构建 撞击 坑 综 合 模型 的 基础 。 


m 


w ox 


析 ， 以 探讨 形态 指标 因子 两 两 之 间 的 量化 关 


重点 。 


理 特 征 的 各 因子 之 间 相 关 关 系 的 密切 程度 ， 主 要 


f 坑 计算 各 项 指标 ， 进 而 进行 指标 


AA 月 表 撞击 坑 形态 指标 单 相关 系数 表 


Table 4 The simple correlation coefficient of morphological indicators of lunar crater 


度 平 均 大 于 退化 撞击 坑 ， 
EDTN, HENDE 
击 坑 中 的 离散 型 撞击 坑 深 径 比 是 退化 离散 型 撞击 坑 的 2 
音 多 ， 退 化 撞击 坑 由 于 坑 底 被 熔岩 填充 或 坑 壁 发 生 块 移 ， 使 得 撞击 坑 的 深度 变化 远 快 于 直径 变化 ， 
因此 ， 新 鲜 撞 击 坑 的 深 径 比 整体 大 于 退化 撞击 坑 的 深 径 比 。 说 明 描 述 撞击 坑 不 同形 态 特征 的 指标 之 
旨 标 间 的 相关 性 


:一 个 复杂 的 渐变 系统 ， 昌 然 其 内 部 各 要 素 之 间 有 概念 描述 和 量化 表达 的 差 
达 ， 选 取 有 效 地 形 因 子 指标 进行 关联 性 分 析 


"aX 
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因子 的 单 相 、 偏 相关 分 
系 。 各 项 形态 指标 单 相关 系数 计算 结果 如 表 4 所 示 。 


直径 深度 体积 DES KSE felt 坑 壁 坡度 坑 层 起 伏 度 

直径 1 0. 474 0.510 -0.256 -0.209 -0.535 -0. 167 0. 415 
深度 0. 474 1 0.497 -0, 212 0. 014 0. 025 0. 260 0. 723 
体积 0. 510 0. 497 1 -0.266 -0.094 -0.259 -0. 135 0. 431 

ial HE -0. 256 -0. 212 -0. 266 1 0. 180 0. 328 0. 270 -0. 305 
体态 比 -0. 209 0. 014 -0.094 0.180 1 0. 317 0.311 -0.179 
深 径 比 -0. 535 0. 025 -0.259 0.328 0. 317 1 0. 666 -0. 108 
坑 壁 坡度 -0. 167 0. 260 -0.135 0.270 0.311 0. 666 1 -0. 080 

DURERE 0. 415 0. 723 0.431 -0.305 -0.179 -0. 108 -0. 080 1 


YE: 样本 为 1407 个 ， 相 关 性 显著 程度 在 0.01 置信 度 水 平 下 。 


Note: the number of sample is 1407. Significant correlation at a confidence level of 0.01. 


通过 单 因子 相关 性 分 析 结 果 可 知 ， 直 径 与 深 径 比 有 着 较 好 的 相关 性 ， 深 度 与 坑 层 起 伏 度 之 间 的 


相关 性 


对 选择 的 


次 试验 发 现 做 4 次 


程度 的 联系 ， 


较 好 ， 深 径 比 与 坑 壁 坡度 间 也 存在 较 好 的 相关 性 ， 这 说明 撞击 坑 形 


8 个 地 形 因 子 指标 进行 主 成 分 
成 分 分 析 是 合理 的 ， 


此 种 联系 在 多 因子 之 间 表 现 得 


更 为 复杂 ， 有 必要 进 一 
分 析 ， 从 而 确定 最 能 表达 撞 了 
计算 结果 如 表 5、 表 6 所 示 。 
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RI 月 表 撞 击 坑 形 态 指标 主 成 分 分 析 结 


Table5 Principal component analysis results of morphological indicators of lunar crater 


STR ban -Z [i PG EE — XE 


因子 综合 相关 性 分 析 。 
因子 及 组 合 


通过 多 


主 成 分 
HER 
1 2 3 4 
直径 0.784 0.049 0.289 0.288 
深度 0.613 0.705 -0.059 0.025 
体积 0.706 0.237 0.273 0.055 
i E -0.559 0.146 0.204 0.733 
体态 比 -0.386 0.404 0.693 -0.393 
深 径 比 -0.620 0.630 -0.273 -0.074 
坑 壁 坡度 -0.404 0.767 -0.072 0.095 
DURERE 0.707 0.425 -0.283 -0.044 
从 表 5 可 以 看 出 ， 第 一 主 成 分 在 接 击 坑 的 体积 、 直 径 、 坑 唇 起 伏 度 和 深度 上 都 具有 较 大 的 正 载 
荷 ， 而 只 有 深 径 比 在 第 一 主 成 分 内 有 较 好 的 负载 荷 ， 以 上 五 个 指标 从 撞击 坑 的 规模 大 小 及 其 发 育 程 
度 的 视角 来 表征 撞击 坑 的 形态 特征 ， 因 此 ， 第 一 主 成 分 可 以 看 做 是 撞击 坑 规 模 及 其 发 育 的 代表 。 第 
二 主 成 分 与 撞击 坑 坑 壁 平均 坡度 和 深 径 比 省 有 较 大 的 正 相 关 ， 这 两 个 指标 与 撞击 坑 的 坑 层 或 坑 壁 发 
育 密切 相关 ， 这 说 明 第 二 成 分 是 从 撞击 坑 坑 唇 和 坑 壁 发 育 状 况 视 角 来 表达 撞击 坑 形 态 特 征 。 虽 然 在 
第 三 主 成 分 中 体态 比 的 相关 性 得 分 较 低 , 但 体态 比 在 表征 撞击 坑 坑 层 形态 的 狭窄 程度 上 有 独特 优势 ， 
因此 ， 在 描述 撞击 坑 形 态 特 征 中 ， 体 态 比 因子 不 可 或 缺 。 在 第 四 主 成 分 中 ， 圆 度 因 子 上 有 具有 较 大 的 
载荷 ， 而 圆 度 从 撞击 坑 的 坑 层 复杂 程度 这 一 独特 视角 表征 着 撞击 坑 的 形态 特征 ， 因 此 对 撞击 坑 形 态 
特征 的 描述 ， 圆 度 是 具有 独特 优势 的 指标 因子 。 
X6 月 表 撞 击 坑 形 态 指 标 主 成 分 分 析 的 总 方差 解释 表 
Table 6 Total variance interpretation table for principal component analysis of lunar crater 
morphological indicators 
初始 特征 值 提取 平方 和 载 入 
成 分 合计 方差 的 % 累计 %% 合计 方差 的 % 累计 % 
1 2.997 37.466 37.466 2.997 37.466 37.466 
2 1.906 23.824 61.290 1.906 23.824 61.290 
3 0.843 10.542 71.833 0.843 10.542 71.833 
4 0.795 9.937 81.770 0.795 9.937 81.770 
5 0.560 6.994 88.764 
6 0.527 6.584 95.348 
7 0.216 2.697 98.045 
8 0.156 1.955 100.000 
同时 ， 从 表 6 中 不 难 发 现 ， 选 择 4 个 主 成 分 ， 其 累计 页 献 紊 达 81. 77%， 保 障 了 对 后 面 撞击 坑 形 
态 特征 发 育 的 指标 线性 模拟 。 而 从 形态 指标 主 成 分 分 析 的 总 方差 解释 表 中 也 可 看 出 ， 第 三 和 第 四 主 


成 分 的 贡献 率 虽 不 及 第 
可 忽略 的 。 在 第 一 主 成 分 中 ， 影 响 撞 避 
表征 的 二 维 指标 因子 ， 其 
较 多 的 联系 ， 并 非 单一 化 的 关联 关系 。 本 文 基于 多 因子 关联 性 分 析 ， 由 在 寻求 撞击 坑 形 态 发 育 与 由 


他 形态 指标 之 间 的 量化 关系 ， 考 虑 到 多 因子 综合 分 析 的 复杂 性 ， 可 以 将 深 径 比 作为 一 个 主 变量 ， 通 


SET d 


一 和 第 二 主 成 分 ， 但 第 三 主 成 分 和 第 四 主 成 分 对 于 撞击 坑 的 形态 的 影响 是 不 
二 坑 规 模 的 指标 因子 众多 ， 而 在 学 界 ， 深 径 比 作为 撞击 坑 发 育 
击 坑 相对 年 龄 的 判断 有 


民 好 的 指示 作用 一 。 它 与 众多 因子 之 间 有 着 
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方式 来 不 断 增加 控 


剖 因 变量 


丸子 间 关 联 性 分 析 。 各 级 偏 相关 结果 如 表 7 所 示 。 


A 7 月 表 撞击 坑 形态 指标 的 一 级 偏 相 关系 数 表 


Table7 First-order partial correlation coefficient table of morphological indicators of lunar crater 


控制 变量 深 径 比 直径 体积 圆 度 体态 比 坑 壁 坡度 
直径 深 径 比 1.000 0.019 0.234 0.249 0.693 
体积 0.019 1.000 -0.164 0.016 -0.059 
E RE 0.234 -0.164 1.000 0.134 0.239 
体态 比 0.249 0.016 0.134 1.000 0.287 
亢 壁 坡度 0.693 -0.059 0.239 0.287 1.000 
体积 深 径 比 1.000 -0.485 0.278 0.305 0.659 
直径 -0.485 1.000 -0.145 -0.189 -0.115 
Ed RE 0.278 -0.145 1.000 0.162 0.245 
体态 比 0.305 -0.189 0.162 1.000 0.303 
议 壁 坡度 0.659 -0.115 0.245 0.303 1.000 
RRE 深 径 比 1.000 -0.494 -0.188 0.278 0.635 
直径 -0.494 1.000 0.474 -0.172 -0.105 
体积 -0.188 0.474 1.000 -0.048 -0.068 
体态 比 0.278 -0.172 -0.048 1.000 0.277 
议 壁 坡度 0.635 -0.105 -0.068 0.277 1.000 
体态 比 REL 1.000 -0.505 -0.243 0.290 0.629 
直径 -0.505 1.000 0.504 -0.227 -0.109 
体积 -0.243 0.504 1.000 -0.255 -0.112 
到 度 0.290 -0.227 -0.255 1.000 0.229 
DEEDEE 0.629 -0.109 -0.112 0.229 1.000 
坑 壁 坡 | 深 径 比 1.000 -0.576 -0.228 0.206 0.155 
度 直径 -0.576 1.000 0.499 -0.222 -0.168 
本 积 -0.228 0.499 1.000 -0.241 -0.055 
E ju 0.206 -0.222 -0.241 1.000 0.105 
体态 比 0.15 -0.168 -0.055 0.105 1.000 


级 1 
分 别 将 深 径 比 因子 与 其 他 除 


的 影响 。 


me. AALS 


dz 


身 以 外 的 指标 间 关 


因子 不 同 程度 地 对 深 径 


前 相 关 结 果 中 ， 作 为 控制 变量 的 直径 、 体 积 、 圆 度 、 体 态 比 以 及 坑 壁 坡度 ， 


比 与 坑 


联 关系 得 到 进 步 简 化 ， 主要 关系 为 : 直径 、 体积 、 
壁 坡度 两 个 因子 之 间 的 关联 关系 产生 了 相对 较 大 
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同时 ， 通 过 表 7 侦 相 关 分 析 结 果 ， 可 以 得 出 ， 深 径 比 对 于 表征 撞击 坑 的 退化 都 是 较 好 的 评判 指 
标 "9， 而 深 径 比 与 坑 壁 坡度 间 的 关联 性 同时 受 直径 、 体 积 、 圆 度 和 体态 比 的 综合 影响 ， 并 且 偏 相关 
性 较 好 ， 这 说 明 深 径 比 因子 与 其 他 五 个 指标 有 着 密切 的 联系 ， 并 且 可 能 存在 有 综合 的 量化 关系 。 


2.4 复 相关 分 析 


从 计量 地 理学 上 相关 性 分 析 的 意义 上 看 ， 指 标 因子 间 的 单 相 关 分 析 和 侦 相 关 分 析 都 是 探讨 两 两 
刀子 之 间 的 关联 性 。 事 实 上 ， 一 个 要 素 的 变化 往往 受 多 种 要 素 的 综合 作用 和 影响 。 从 表 4 分 析 结 果 
中 可 以 得 知 , 深 径 比 指标 和 深度 指标 都 可 看 作 是 量化 撞击 坑 退 化 研究 中 的 因 变 量 。 因 此 本 文 基于 1407 
个 撞击 坑 样 本 数据 ， 探 讨 最 终 因 变 量 的 选择 以 及 因 变 量 与 其 他 因子 之 间 的 复 相关 关系 ， 以 量化 各 指 
标 因子 对 于 因 变 量 的 综合 影响 。 复 相关 回归 统计 结果 如 表 8 所 示 。 

表 8 月 表 撞 击 坑 形 态 指标 的 回归 统计 结果 表 


Table8 The regression statistics table of morphological indicators of lunar crater 


F 复 相关 系数 R RV 标准 误差 观测 值 
~ 深度 435.513 0.780 0.609 890.857 1407 
深 径 比 492.037 0..798 0.637 0.0392 1407 


结果 可 得 ， 深 径 比 因子 与 其 他 5 个 因子 之 间 的 相关 性 更 好 ， 且 高 度 显著 ， 因 此 选择 深 径 比 作为 
探索 各 指标 关联 性 的 因 变量 。 而 深 径 比 与 直径 、 体 积 、 圆 度 、 体 态 比 和 坑 壁 坡度 因子 之 间 的 相关 性 
也 进一步 证 明 各 指标 作为 撞击 坑 形 态 特 征 的 指示 因子 ， 存 在 着 密切 的 数学 联系 ， 可 进行 多 变量 的 综 
合 回 归 模型 分 析 。 


2.5 基于 多 因子 的 量化 关系 模型 构建 


虽然 指标 间 的 关联 性 分 析 已 经 说 明 深 径 比 与 直径 、 体 积 、 圆 度 和 体态 比 、 坑 壁 坡 度 5 个 指标 因 
c 子 存在 着 综合 的 数学 关系 。 但 从 地 理 建 模 角 度 来 看 ， 各 自 变 量 应 当 对 于 因 变 量具 有 很 好 的 可 解释 性 
r- 和 数学 依赖 性 ， 并 且 自 变量 之 间 应 当 保持 一 定 的 独立 性 。 
Q 深 径 比 因子 与 直径 、 坑 壁 坡度 因子 之 间 的 相关 性 较 好 ， 而 圆 度 因子 、 体 态 比 因子 与 其 他 指标 因 
子 之 间 的 相关 性 中 ， 与 深 径 比 的 关联 性 最 好 。 并 且 除 深 径 比 之 外 的 其 他 因子 从 不 同 视角 对 于 撞击 坑 
形态 都 具有 较 好 的 地 学 表征 意义 。 从 单 相关 分 析 结果 中 看 ， 在 作为 五 个 自 变量 因子 之 间 的 相关 性 比 
较 一 般 ， 其 因子 之 间 的 独立 性 较 好 。 

综 上 所 述 ， 本 文 将 以 深 径 比 为 因 变量 CY), 以 直径 Q), AR OG) 、 圆 度 Oe) 、 体 态 比 XO 
以 及 坑 壁 坡度 (Xs) 为 自 变量 ， 构 建 多 元 线性 回归 的 量化 关系 模型 。 采 用 逐步 (SWR) 回归 模拟 发 
现 ， 各 变量 之 间 存在 很 强 的 线性 函数 关系 ， 其 方程 为 公式 OD : 


Y = (7.941E - 007) x x, + (8.063E - 016) x x , +0.082xx, +0.012xx, +0.007 xx, -0.042 (1) 


2.6 模型 的 统计 检验 


本 文 从 模型 的 拟 合 优 度 、 方 程 的 显著 性 、 变 量 显著 性 与 参数 的 置信 区 间 四 个 方面 对 模型 进行 统 
计 检 验 。 拟 合 模 型 的 复 相关 系数 为 0.798， 回 归 模 型 的 概率 为 0.00， 明 显 小 于 显著 性 水 平 0.05。 模 
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型 的 拟 合 优 度 达 0.637， 调 整 后 的 拟 合 优 度 达 0. 639。 回 归 模 型 的 方差 分 析 结 果 见 表 9。 经 过 显著 性 
检验 ， 在 0. 05 的 置信 水 平 下 查 F 分 布 表 可 知 : F k,n-k-UDJ=EF,.(51400D)= 2.21， 在 回归 拟 合 结果 中 ， 
F=492.037>>Fao KIJE, Æ 0. 05 的 置信 水 平 下 ， 回 归 模 拟 函 数 是 极 显 车 的。 


表 9 撞击 坑 多 元 线性 回归 模型 的 方差 分 析 


Table 9 Analysis of variance of multiple linear regression model of crater 


平方 和 df 均 方 F Sig.F 
可 归 3.787 5 0.757 492.037 0.00^ 
残 差 2.156 1401 0.002 
总 计 5.943 1407 


3 小 结 与 展望 


本 文 从 撞击 坑 的 尺寸 特征 、 发 育 特 征 等 角度 总 结 了 前 人 对 撞击 坑 形 态 特征 的 基本 认识 和 量化 研 
究 成 果 ; 利用 100 米 分 辩 率 的 DEM 数据 ， 通 过 对 撞击 坑 各 形态 因子 分 析 筛 选 给 出 了 撞击 坑 形 态 的 多 
元 线性 回归 拟 合 模型 ， 这 种 拟 合 关系 反映 了 撞击 坑 形 态 特征 的 指标 因子 传递 出 的 撞击 坑 发 育 的 数学 
机 理 ， 揭 示 了 撞击 坑 形 态 发 育 指标 因子 之 间 的 内 在 联系 ， 而 不 再 是 单一 的 因子 量化 特征 表达 。 通 过 
对 描述 月 表 撞 击 坑 形 态 参数 的 指标 的 关联 性 分 析 和 量化 模拟 研究 ， 可 获得 以 下 结论 : 

(1) 在 文 [13] 作 者 对 撞击 坑 分 类 的 基础 上 ， 将 撞击 坑 划 为 6 种 类 型 ,对 选取 的 1407 个 撞击 坑 进 
行 量 算 ， 新 鲜 撞 击 坑 的 深 径 比 平均 值 为 0.187， 退 化 撞击 坑 的 深 径 比 平均 值 为 0.0771， 其 他 撞击 坑 形 
态 参 数 指标 的 差异 如 表 2 示 。 

(20 对 8 个 撞击 坑 形 态 指标 进行 偏 相关 和 复 相 关 的 分 析 , 确定 以 深 径 比 为 因 变 量 , 直径 、 体 积 、 
圆 度 、 体 态 比 和 坑 壁 坡度 为 自 变量 构建 指标 因子 的 量化 模拟 ， 得 到 其 线性 拟 合 公式 〈 公 式 1) Bi 
型 的 拟 合 优 度 达 0. 637, 调整 后 的 拟 合 优 度 达 0. 639， 说 明 本 文 的 拟 合 模 型 能 够 较 好 地 表达 撞击 坑 形 
态 发 育 特征 。 

(3) 模型 的 复 相关 系数 为 0.798， 在 0. 05 置信 水 平 的 显著 性 F=492. 037>>F。， 说明 构 建 的 模型 
显著 性 较 好 。 可 用 于 月 表 接 击 坑 形 态 发 育 特征 的 预测 和 估计 。 

此 外 ， 本 文 还 对 所 选择 的 指标 因子 进行 非 线 性 拟 合 ， 深 径 比 与 直径 、 坑 壁 坡度 的 拟 合 优 度 分 别 
为 0. 4715、0. 4524， 而 深 径 比 与 体积 、 圆 度 、 体 态 比 因子 的 拟 合 优 度 都 小 于 0. 15， 也 说 明 深 径 比 与 
其 他 指标 因子 之 间 不 适宜 多 元 非 线性 拟 合 。 

本 文 基于 撞击 坑 形 态 特 征 ， 实 现 多 因子 量化 关系 建 模 来 表达 撞击 坑 形 貌 发 育 的 演化 机 理 ， 揭 示 
描述 撞击 坑 形 貌 的 地 形 因 子 信息 与 月 球 空间 环境 系统 之 间 规 律 的 过 程 ， 对 于 多 因子 联合 定量 表达 的 
撞击 坑 形 貌 信息 建 模 研 究 和 撞击 坑 形 貌 形 态 发 育 研 究 有 重要 的 实践 意义 。 

然而 ， 在 利用 建 模 的 理论 和 方法 实现 撞击 坑 形 态 数字 地 形 信息 挖 气 的 方面 ， 还 需 开 展 大 量 的 工 
作 。 在 指标 选取 上 要 进行 更 多 的 尝试 ， 以 增强 模型 构建 的 稳定 性 。 提 取 全 月 撞击 坑 形 态 指标 数据 ， 
对 模型 进行 全 月 撞击 坑 数 据 的 验证 ,探讨 以 模型 为 依据 , 结合 月 表 岩 性 分 类 研究 ! 的 撞击 坑 类 型 分 类 
研究 是 未 来 的 研究 工作 。 
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Abstract: Based on the surface DEM data of 100m resolution obtained by the laser altimetry method, 


1407 different types of impact craters are Selected to extract the topographic factor indicators of the 
crater morphology from multiple perspectives, and the morphological characteristics of the crater are 
analyzed. On the basis of the method of mathematical statistics and geographic modeling, the 
relationship between topographic factors is mainly discussed, and a multi-factor based quantitative 
model is further constructed to explain the intrinsic relationship between topographic factor 
information and crater morphology development. The results show that the depth-diameter ratio (Y) 
and diameter (X4), volume (X3), circularity (X5), posture ratio (X4) and crater wall slope (X5) can better 
characterize the morphological characteristics of the crater. Multiple linear regression fitting was used 
to construct the functional relationship between the indicators, and the significance of the fitted model 
at the 0.05 confidence level was F=492.037>>F. The significance of the fitted model was excellent, 
and its simulation function relationship was fitted. Goodness is 0.637. The simulation method is of 
great significance in the study of the morphological development of the lunar impact crater. It is a 
useful practice of digital terrain simulation and information mining of the crater feature. 


Key words: crater morphology; Factors of terrain; Crater development; Quantitative simulation 


